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Continental ice sheets typically sculpt landscapes via erosion; under certain 
conditions, ancient landscapes can be preserved beneath ice and can survive 
extensive and repeated glaciation. We used concentrations of atmospherically 
produced cosmogenic beryllium-10, carbon, and nitrogen to show that ancient soil 
has been preserved in basal ice for millions of years at the center of the ice sheet at 
Summit, Greenland. This finding suggests ice sheet stability through the 
Pleistocene (i.e., the past 2.7 million years). The preservation of this soil implies that 
the ice has been non-erosive and frozen to the bed for much of that time, that there 
was no substantial exposure of central Greenland once the ice sheet became fully 
established, and that preglacial landscapes can remain preserved for long periods 
under continental ice sheets. 


A diverse set of geochemical records has been developed from ice re- 
covered in the 3054-m Greenland Ice Sheet Project 2 (GISP2) core. 
These data provide a detailed history of climate and ice dynamics 
stretching back over 100,000 years. The lowermost 13 m of the core, 
which contain high concentrations of sediment, are less well studied. 
Here, we used measurements of atmospherically produced (meteoric) 
^^Be, carbon, and nitrogen in sediment extracted from this silty ice zone 
to better understand the source and age of this material and to infer the 
erosional and thermal history under the center of the Greenland Ice Sheet 
(GIS) since its inception millions of years ago. 

Ice likely was present on Greenland starting ~38 million years ago 
(Ma), expanding after 14 Ma with glaciations becoming more frequent 
after ~7 Ma and full coverage of Greenland by ice occurring at ~2.7 Ma 
(7, 2). Fossils in north Greenland suggest a partially forested tundra 
landscape some time before the growth of the current ice sheet (3). The 
more recent history of the GIS (the past 100,000 years) is constrained by 
ice cores as well as geologic mapping (4), cosmogenic nuclides (5), geo- 
chemical proxies ( 6 ), sea level records (7), and dating of sediment 
(8). 

Glacial erosion rates are controlled by thermal conditions at the bed. 
Ice colder than the pressure melting point is frozen to the bed (cold- 
based) and non-erosive. Ice warmer than the pressure melting point 
(warm-based) can erode the bed. Using heat-flux data, one model sug- 
gests that >75% of the present GIS is warm-based and thus erosive (9); 
however, the present-day thermal state of the GIS is known with certain- 
ty only where drilling exposes the bed. 


The rate at which the GIS erodes its 
bed and the spatial distribution of sub- 
ice erosion are uncertain. North Atlan- 
tic sediment volume estimates suggest 
that, on average, -100 m of rock were 
removed from glaciated areas of North 
America and Greenland after sustained 
glaciation began at 2.7 Ma ( 10 ). This 
corresponds to an average erosion rate 
of -40 m per million years — a value 
that masks variability over space and 
time. 

All GIS cores that have been drilled 
to or near to the bed contain silty basal 
ice (77). This ice can contain organic 
material (72), microorganisms, and 
high gas concentrations (73) and can 
display a stable water isotope signature 
( 14 ) interpreted as indicating warmer 
conditions when the basal ice formed 
(77), perhaps in the absence of an ice 
sheet (75). 

At the bottom of the GISP2 ice core 
(Summit, Greenland; Fig. 1), 13 m of 
silty ice [average 0.4% solids ( 16 )] 
overlies 48 cm of diamict containing 
erratic lithologies and underlain by 
granitic rock (77, 16 , 17 ). The basal ice 
has not been dated directly, but the 
composition of dissolved gases sug- 
gests an age of at least 237,000 years 
( 18 ). Ice at Summit is currently frozen 
to the bed (-9°C, 6°C below the pres- 
sure melting point) (79), and models 
suggest that it has remained frozen to 
the bed for at least the last several gla- 
cial cycles ( 20 ). 

To understand the origin and history of the sediment in basal GISP2 
ice, we measured the concentration of meteoric ^^Be adhered to sand and 
silt extracted from 17 samples of GISP2 silty ice. This isotope is contin- 
uously produced by cosmic rays in the atmosphere, is delivered to 
Earth’s surface by precipitation and dry fall, and adheres to sediment. It 
has been used to date soils and infer erosion rates (27). We also meas- 
ured organic carbon (OC) and total nitrogen (TN) in these same samples. 
To provide an analog for the GIS data, we measured a meteoric ^^Be 
profile in a permafrost tundra soil developed on a mid-Pleistocene gla- 
cial surface in Alaska (22). Methodological details and data tables are 
provided in the supplementary materials. 

All silty ice sediment samples contained high concentrations of me- 
teoric ^^Be and detectable levels of OC and TN (Fig. 2). Sediment in an 
11 -cm sample of ice just above the bed had the highest concentration of 
meteoric ^^Be (3.8 x 10^ atoms g~^). Other sediment samples contained 
0.61 X 10^ to 1.02 X 10^ atoms meteoric ^^Be g“\ In general, we found 
that ^^Be and OC concentrations decrease with sample distance above 
the bed. This decrease is consistent with the mixing of two components: 
(i) locally eroded soil material with both high meteoric ^^Be and OC 
concentrations, and (ii) eroded rock or sediment containing little meteor- 
ic ^^Be or OC. The second component could be rock flour derived from 
outcrop abrasion or till sourced from below the depth to which meteoric 
^^Be has penetrated soil, regolith, and rock. 

The geochemical data suggest that sediment incorporated into the silty 
ice zone was derived from erosion of a stable cold-region landscape. 
^^Be concentrations measured in silty ice sediment are similar to those 
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measured in long-lived soils (27). Compared to 541 ^^Be concentrations 
measured in soil (27), the maximum ^^Be concentration measured in 
GISP2 sediment ranks in the 60th percentile (Fig. 3 A). Sediment ex- 
tracted from GISP2 silty ice contains 0.3 to 1.7% OC and 0.03 to 0.14% 
TN. The OC/TN regression indicates a C/N ratio of ~10, consistent with 
values reported for mineral-rich subsoils from boreal regions and cry- 
oturbated tundra (23, 24). Both OC and TN concentrations are positively 
correlated with ^^Be concentration and all decrease above the bed, indi- 
cating that ^^Be, C, and N are derived from the same source, most likely 
a former tundra soil. 

We can estimate the inventory of meteoric ^^Be in the soil before 
burial beneath the GIS. Such a calculation is viable because soil general- 
ly retains ^^Be and because ^^Be soil inventories — at least in temperate 
regions where sufficient data exist to make this comparison — are well 
correlated on a global scale with the maximum measured concentration 
of ^^Be in the soil profile (Fig. 3B) (27). Taking the maximum measured 
meteoric ^^Be concentration (3.8 x 10^ atoms g^^) at face value suggests 
an inventory of 6.7 x 10^^ atoms meteoric ^^Be cm~^ (Fig. 3B). However, 
if we presume that the soil from which the silty ice was derived predates 
the GIS and is thus older than 2.7 million years (~2 half-lives of ^^Be), 
then the concentration when the soil was first buried by ice, after correct- 
ing for radioactive decay of ^^Be, was higher by a factor of more than 4 
(-1.5 X 10^ atoms meteoric ^^Be g^^). The corresponding soil inventory 
when the GIS developed was -3.8 x 10^^ atoms meteoric ^^Be cm“^. 
Such concentrations and inventories, while high, are reasonable. The 
tundra soil from the North Slope of Alaska, developed on a till plain that 
has been stable for between 150,000 and 750,000 years (22), contains up 
to 1.1 X 10^ atoms meteoric ^^Be g~^ and has an inventory of 4.0 x 10^^ 
atoms meteoric ^^Be cm~^, consistent with the meteoric ^^Be concentra- 
tions in the silty ice zone and the inferred meteoric ^^Be inventory (Fig. 
4). 

The duration over which this ancient soil developed in Greenland be- 
fore expansion of the GIS can be roughly estimated using the delivery 
rate of meteoric ^^Be to the GISP2 site [3.5 x 10^ atoms cm~^ year^^ in 
the upper, silt-free ice of the GISP 2 core (25)]. With no erosion and 
assuming the Holocene meteoric ^^Be deposition rate, the soil material in 
the silty ice zone was subaerially exposed for at least 200,000 years (no 
decay correction) and perhaps more than 1 million years (decay- 
corrected). Even using a deposition rate several times higher, assuming a 
weakened polar vortex and increased precipitation during past warm 
periods, the soil incorporated in GISP2 basal silty ice requires far longer 
exposure (on the order of 10^ years) than can be accounted for by a few 
millennium-long periods of exposure during extreme interglacials [such 
as marine isotope stage (MIS) II]. Therefore, the data are most con- 
sistent with soil formation prior to the existence of the present GIS. On 
the basis of similarly high meteoric ^^Be concentrations, ancient pregla- 
cial regolith has been identified on uplands and in glacial sediments 
preserved on Baffin Island in arctic Canada (26). 

Other means of generating high meteoric ^^Be concentrations in sed- 
iment extracted from the GISP2 silty ice zone do not appear plausible. It 
is unlikely that the ^^Be we measured was concentrated from basal 
meltwater. Ice in the GISP2 core contains on average 1.4 x 10"^ atoms 
meteoric ^^Be g~^ (25), four orders of magnitude less than measured in 
the silty ice. Stable water isotope data from basal ice of the nearby 
Greenland Ice Core Project (GRIP) core (79) fall on the meteoric water 
line and are inconsistent with large-scale melting that would be needed 
to concentrate ^^Be. The positive, linear correlation between OC and 
^^Be (Fig. 2) is also inconsistent with delivery of meteoric ^^Be from 
meltwater. Long and/or warm interglacials (MIS 5e, 9, and II) are un- 
likely to have melted the entire GIS, and even if they did, the peak of 
each interglacial is so short (thousands of years) that insufficient meteor- 
ic ^^Be would be delivered to any exposed soil at Summit to build up 
even a few percent of the meteoric ^^Be inventory we infer. On the basis 


of the marine isotope record (27), other interglacials were less intense 
than MIS 5e, 9, and 11 and thus even less likely to melt the entire GIS. 
In situ production of ^^Be by the interaction of cosmic rays (specifically 
muons) with sediment under the ice sheet, which is several kilometers 
thick, generates an inventory of only a few atoms ^^Be g^^ over a period 
of 2.7 million years (28). 

Previous but imprecise measurements of in situ ^^Be, ^^Al, and ^^Cl 
in rock collected below the silty ice zone of the GISP2 core are con- 
sistent with short exposure of that rock (at most a few thousand years) 
around 0.5 + 0.2 Ma, possibly during the long MIS II interglaciation 
(29). If there were a short period of exposure during MIS 11, then the 
pre-Pleistocene soil below GISP2 was not fully eroded during that expo- 
sure. Any in situ produced ^^Be, ^^Al, and ^^Cl accumulated in rock dur- 
ing brief exposure at MIS 1 1 could have been produced under a cover of 
uneroded, pre-Pleistocene soil, or the soil we sampled could have been 
transported from nearby. The ^^Be measurements presented here demand 
a much longer period of subaerial exposure than the data of Nishiizumi 
et al. (29), many tens to a few hundreds of thousands of years, a duration 
of exposure only possible at Summit before the expansion of the GIS at 
-2.7 Ma. 

The continued presence of soil in GISP2 ice for several million years 
after formation of the GIS indicates extremely low rates of sub-ice ero- 
sion and slow rates of horizontal ice advection away from Summit. Such 
low shear rates are consistent with both long-term stability of the GIS ice 
divide location and a frozen bed for most, if not all, of the past 2.7 mil- 
lion years. The existence of ancient soil implies low sub-ice erosion 
rates; however, the 13-m thickness of the silty ice zone suggests erosion, 
entrainment, and mixing of the preglacial soil. Although persistence of 
this ancient soil precludes substantial erosion at the bed (16), it is possi- 
ble that cold-based ice incorporated some basal debris (30). The flux of 
silty ice past the GISP2 site is likely dominated by deformation within 
the ice rather than by sliding (77). If warm-based ice was present at 
Summit for short periods of time in the past, it did little more than mix 
sediment from the bed into the overlying ice. 

Meteoric ^^Be data constrain both the preglacial landscape history 
and the Quaternary history of Summit, Greenland. ^^Be and OC data are 
consistent with inferences, based on the stable oxygen isotopic composi- 
tion of ice and the concentration and isotopic composition of gases in the 
silty ice zone, that the GIS overran a local snowfield or ice that itself 
overlay a soil (75). Since the GIS formed, the soil has been preserved 
and only slowly eroded, implying that an ancient landscape underlies 
3000 m of ice at Summit. These new data are most consistent with con- 
tinuous cover of Summit by ice for the entire Quaternary, with at most 
brief exposure and minimal surface erosion during the warmest or long- 
est interglacials. 
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Fig. 1. Photograph of GISP2 ice core silty ice zone; the inset shows core stratigraphy ( 15 ). Thick gray lines indicate intervals 
of the ice core that were analyzed for meteoric ^°Be, OC, and TN. Core depth shown in meters by labels. Image from NOAA. 
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Fig. 2. Geochemical and isotopic data for GISP2 basal silty ice. (A) Meteoric ^°Be concentration of silt with core depth; 
horizontal axis has gap and variable scaling to allow plotting of the high-concentration lowermost sample. (B) Water oxygen 
isotopes in ice ( 14 ). (C) Organic carbon (OC) measured in acidified silt. (D) Organic carbon/total nitrogen ratio measured in 
acidified silt. (E) Positive linear relationship between organic carbon and meteoric ^°Be in silt samples. Data in table S1 . 
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Fig. 3. GISP2 silt contains high concentrations of meteoric ^°Be. (A) Histogram of 541 measurements of meteoric ^°Be made 
in soils (21) with maximum meteoric ^°Be measured in GISP2 silt and decay-corrected concentration assuming burial for two 
half-lives of ^°Be (~2.7 million years). (B) Maximum meteoric ^°Be concentrations in soil profiles are well correlated to soil 
meteoric ^°Be inventories (21). Using measured and decay-corrected maximum meteoric ^°Be concentrations in GISP2 silty 
ice suggests that meteoric ^°Be inventory in pre-GIS soil was at least 6.7 x 10^° atoms cm“^ and likely 3.8 x 10^^ atoms cm"^. 
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Fig. 4. Depth profile of meteoric ^°Be measured in a stable 
soil profile developed on glacial sediment deposited 
>150,000 years ago, North Slope, Alaska. Prominent 
subsurface bulge is indicative of a stable soil profile (21). 
Maximum measured and decay-corrected concentrations of 
meteoric ^°Be measured in GISP2 silt are overlain. Data in 
table S2. 
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